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ABSTRACT 



Context. Intermediate-mass black holes (IMBHs) fill the gap between stellar-mass black holes and supermassive black holes 
(SMBHs). The existence of the latter is widely accepted, but there are only few detections of intermediate-mass black holes (10^ 
- 10^ Mg) so far Simulations have shown that intermediate-mass black holes may form in dense star clusters, and therefore may still 
be present in these smaller stellar systems. Also, extrapolating the M. - a scaling relation to lower masses predicts intermediate-mass 
black holes in systems with a ~ 10-20 km/s such as globular clusters. 

Aims. We investigate the Galactic globular cluster NGC 5286 for indications of a central intermediate-mass black hole using spectro- 
scopic data from VLT/FLAMES*, velocity measurements from the Rutgers Fabry Perot (RFP) at CTIO, and photometric data from 
HST/ACS. 

/Wefhods. We compute the photometric center, a surface brightness profile, and a velocity-dispersion profile. We run analytic spherical 
and axisymmetric Jeans models with different central black-hole masses, anisotropy, mass-to-light ratio, and inclination. Further, we 
compare the data to a grid of A'-body simulations without tidal field. Additionally, we use one A'-body simulation to check the results 
of the spherical Jeans models for the total cluster mass. 

Results. Both the Jeans models and the A'-body simulations favor the presence of a central black hole in NGC 5286 and our detec- 
tion is at the 1- to 1.5-cr level. From the spherical Jeans models we obtain a best fit with black-hole mass M. = (1.5 ±1.0) • 10^ Mq. 
The error is the 68% confidence limit from Monte Carlo simulations. Axisymmetric models give a consistent result. The best fitting 
A'-body model is found with a black hole of 0.9% of the total cluster mass (4.38 ± 0.18) • lO'' Mq, which results in an IMBH mass of 
M. = (3.9 ± 2.0) • 10' Mq. Jeans models give values for the total cluster mass that are lower by up to 34% due to a lower value of M/L. 
Our test of the Jeans models with A'-body simulation data shows that the discrepancy in the total cluster mass has two reasons: The 
influence of a radially varying M/L profile, and underestimation of the velocity dispersion as the measurements are limited to bright 
stars, which have lower velocities than fainter stars. We conclude that detection of IMBHs in Galactic globular clusters remains a 
challenging task unless their mass fractions are above those found for SMBHs in nearby galaxies. 

Key words. Black hole physics - globular clusters: individual: (NGC 5286) ~ stars: kinematics and dynamics 



1. Introduction 

The existence of supermassive black holes (SMBHs) is widely 
accepted. Mass estimates come either from kinematic measure- 
ments, or from X-ray and radio luminosities, which assume 
certain accretion rates. However, the formation and growth of 
SMBHs is not well understood, although it is somehow linked 
to the evolution of the environment. The mass of a central black 
hole is correlated to the mass of the bulge M^ of its host system 
(e.g. Haring & Rix|2004)i, and to the stellar velocity dispersion 



cr of the bulge (e.g.'Gebhar dt et al.|2000a[ [Ferrarese & Merritt| 
[2000 , Gultekin et al. 2009 ). Systems with higher Mh and higher 
cr contain more massive central black holes. This suggests a con- 



* Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile 
(085.D-0928) 



nection between the formation and evolution of central black 
holes and their host systems. 

There are several detections of quasars at a redshift higher 
than 6 (e.g. [Fan et al.|200T] [Willott et al.lSoTOl [Mordock et aL 
|2011[ ). The masses for the central black holes as inferred from 
the luminosity are of the order of ~ 10^ Mg. These high masses 
cannot be explained by accretion onto a stellar-mass black hole, 
as not enough time has elapsed at z ~ 6: the accretion rates would 
greatly exceed the Eddington limit. In order to build a black 
hole of 10'' Mq at z ~ 6, the black-h ole see d requires a mass of 
10^-10^ Mq, as shown by Volonteri (2010). Black holes in this 
mass range are called intermediate-mass black holes (IMBHs) 
and must have formed in the young universe. 

There are theories for the formation of IMBHs from 
Population III stars (Madau & Rees [2001 Bromm & Loeb 



2003), gas-dynamical processes (Loeb & Rasio 



stellar-dynamical instabilities in star clusters (Ebisuzaki et al. 



19941, and from 
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|2001| l. The latter theory is of particular interest for globular clus- 
ters. In very compact star clusters the most massive stars sink to 
the core due to mass segregation. If the cluster is very dense, 
this occurs befo re the stars explode an d evolve to compact rem- 



tent with an IMBH of ~ 2 ■ IO'^Mq. Miller- Jones et al. (2012 1 



nants (~ 10^ yr, Ebisuzaki et al. 12001 1. In a dense core the mas 



sive stars undergo runaway stellar mergers. After several col 

lisions, a very massive star (VMS) with > 100 Mq forms and 

grows quickly (Portegies Zwart et al. 1999|l. Yungelson et al. 



confirmed the detection of X-ray emission but detected no ra- 
dio emission in Gl. They concluded that the detected emission 
comes more likely from a low-mass X-ray binary (LMXB) with 
a stellar-mass black hole than from an IMBH. 

Recent studies found no significant radio sources in the cen- 
ter of OjC entauri ( |Lu & Kong][20n) and NGC 6388 ( |Cseh 



( 2008 1 used an evolutionary code to compute the mass loss for 
VMSs of initial solar metallicity (Z = 0.02). They found that the 
mass loss due to stellar winds is so high that it is unlikely that the 
VMS forms an IMBH. Belkus et al. (2007) studied the evolution 
of VMSs at different metallicities, and found that low-metallicity 
stars with Z - 0.001 suffer less mass loss. Their final mass is by 
a factor of ~ 2-3 larger than for stars with a metallicity of Z 
= 0.02. These VMSs can evolve to IMBHs with some 10^ Mq. 
The central black hole continues accreting nearby stars, until the 
reservoir of massive stars is exhausted. 

If a cluster is not dense enough, stellar evolution is faster 
than mass segregation. Massive stars evolve to stellar remnants 
before reaching the center. [Miller & Hamilton| ( |2002 [) proposed 
that successive mergers of stellar remnants and accretion can 
form an IMBH in a globular cluster. Recent models examined 
the retention fraction of stellar black holes in star clusters and 
merger rates due to gravitational radiation and found different 
results. The retention fraction for black holes binaries is only 
1% to 5% acc ording to [Moody & Sigurdsson| ([2009]). [Miller &| 



et al. 2010[l. However, mass estimates from upper limits on radio 



emission require assumptions on the gas properties and accretion 



model. Pepe & Pellizza 



( 2013p showed that the accretion rate in 



globular clusters depends on the mass of the cluster and its veloc- 
ity dispersion. Further, Perna et al. (|2003[) found tha t the widely 

_ .._--. ■" ■ ■ ■ an 



T944) l is 



used Bondi-Hoyle accretion rate (Bondi & Hoyle 

upper limit of the true mass accretion, which can be lower by 

several orders of magnitude. The fraction of the Bondi accretion 

rate at which accretion occurs is highly uncertain ( Maccarone] 

|2004| . 

There are further examples for black holes in extragalac- 
tic globular clusters. [Zepf et al. (2008) measured broad emis- 
sion lines in the globular cluster RZ 2109 in the Virgo elliptical 
galaxy NGC 4472 and conclude that, in combination with the 
measured X-ray luminosity, this can be explained by the pres- 
ence of a ~ IOMq black hole. A globular cluster in NGC 1399 
shows broad emission lines and X-ray emission that can be mod- 
eled with the tidal disruption of a horizontal branch star by a 
massive (50-lOOMo) black hole ( [Clausen et al.|2012 !). 



Davies ( 2012) found that nearly all black holes are kicked out 
during mergers. On the other hand, Morscher et al. (2013|) came 



Baumgardt et al. ( 2005 1 performed A^-body simulations and 



to the result that less than 50% of the black holes are dynamically 
ejected, and the total merger rate is ~ 1 per Gyr and globular 



cluster. Also O'Leary et al. (2006) found significant probability 
of black hole growth. However, runaway mergers of black holes 
to form an IMBH seem to be less likely for globular cluster-sized 
systems. 

Globular clusters orbiting around galaxies will sink to the 
galactic center due to dynamical friction ( Ebisuzaki et al.|200T i. 
While stars in the outskirts of the clusters are stripped away from 
the galaxy's tidal field and settle in the galactic bulge, the central 
IMBHs can be carried to the center of the galaxy. They can form 
binary IMBHs. Some IMBHs are ejected in three-body interac- 
tions, and some IMBHs merge through gravitational radiation 
into more massive black holes, and form a SMBH by successive 
mergers ( Ebisuzaki et al. 2001) . If the M, - cr scaling relation 
holds down to the low-mass end, globular clusters, which have 
velocity dispersions of ~ 10 - 20km/s, are a possible environ- 
ment for intermediate-mass black holes. 

INoyola et al. (2008) found kinematic signatures for a black 
hole of (4.0 ;^j y"") ■ 10"* Mq in wCentauri from the radial veloc- 
ity dispersion and isotropic Jeans models. Their result was chal- 
lenged by 'van der Marel & Anderson (2010 1, who used HST 
proper motion measurements and determined an upper limit 
of 1.8- 10''Mq. The groups used different center coordinates, 
which is probably one of the main reasons for this discrepancy. 
With the kinematic center from |Noyola et al.| ( [2010| ), [Jalali et al.[ 
( 2012| ) performed direct A^-body simulations and found a central 
IMBH with ~ 5 ■ lO'* M©, which is consistent with the spherical 
isotropic models of Noyol a et al.| ( [201^ . Another good candi- 
date with kinematic measurements is NGC 6388. [Liitzgendorf 
et al. ( 2011[ ) used integral field unit data and Jeans models to de- 
termine a black-hole mass of (1.7 + 0.9) ■ 10"^ M©. The most mas- 
sive globular cluster in M31, Gl, shows kinematic signatures 
([Gebhardt et al. 2005) as well as X-ray (Pooley & Rappaport 



[2006] ) and radio emission ( |Ulvestad et al.||2007j ), all consis- 



ascertain that globular clusters with an IMBH have a rather large 
core, and a shallow cusp in the surface brightness profile. The 
velocity-dispersion profile rises towards the center, but this is 
difficult to detect if only the brightest stars of a cluster are ob- 
served. Noyola&Gebhardt (2006) obtained a surface brightness 
profile of NGC 5286 that shows a shallow cusp, and for this rea- 
son we set out to search for kinematic evidence for an IMBH. 

NGC 5286 is in the halo of the Milky Way at a distance of 
about 1 1.7 kpc ([Har ris 1996). Tableflllists some properties of this 
cluster. [Milone et al.| ( |2012| ) found some foreground stars in the 
region of NGC 5286. The fairly bright (my = 4.6) spectroscopic 
binary star M Centauri (HR 5172) is at a distance of about 4' 
from the center of NGC 5286. [Zoroto vie et al.[ ( f2009 ) compared 
the CMD of NGC 5286 with M3 and found that NGC 5286 is 
about (1.7 ±0.9) Gyr older than M3. It may belong to the old- 
est globular clusters in our Galaxy. [Crane et al.[ ( [2003| ) argued 
that NGC 5286 is associated with a ringlike structure near the 
Galactic anticenter, based on its position and kinematics. This 
ringlike filament of stars around the Milky Way, the so-called 
Monoceros ring, can be the tidal stream of a disrupted dwarf 
galaxy, with the nucleus in Canis Major (jForbes et al.|2004i. 
But the existence of the Canis Major dwarf galaxy is still un- 
der debate. However, it is possible that NGC 5286 is originally a 
globular cluster of the Canis Major dwarf galaxy and not of the 
Milky Way. 

This paper is organized as follows: In Sectionl2]we describe 
the photometric analysis, which includes the center determina- 
tion and calculation of a surface brightness profile. This gives 
information on the stellar distribution and density. The observa- 
tion and data reduction of the spectroscopic FLAMES data are 
outlined in Section [3] With these data we calculate a velocity 
map and a velocity-dispersion profile for the central 25", as de- 
scribed in Section]?] For the outer part of the profile we use a data 
set obtained from the Rutgers Fabry Perot, and we investigate it 
for indications of rotation. Sectionl5]and Sectionl6]introduce the 
results of the Jeans models and A^-body simulations, which are 
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Table 1. Properties of the globular clus ter NGC 5286. 
Refere nces: GolO = [Goldsbury eTill ([20T0| , Hal O ^ [Harris 
([19961 2010 edition), Zi80 ^ |Zinn[ (|1980|l, ZolO = [Zorotovic 



Table 2. Center coordinates (J2000) 



^taL( 



20T0I 1, Wh87 ^ [White & Shawl| ( |1987| l 



Parameter 



Value 



Reference 



RA a (J2000) 


13h 46m 26.81s 


GolO 


DEC 6 (J2000) 


-51° 22' 27.3" 


GolO 


Galactic longitude 1 


311.6142° 


GolO 


Galactic latitude b 


+ 10.5678° 


GolO 


Distance from the sun 


11.7 kpc 


HalO 


Distance from the Galactic center 


8.9 kpc 


HalO 


Foreground reddening E(B-V) 


0.24 


Zi80 


Absolute visual magnitude My 


-8.74 


HalO 


Metallicity [Fe/H] 


-1.67 


ZolO 


Projected ellipticity e = l-(b/a) 


0.12 


Wh87 


Heliocentric radial velocity v, 


(57.4 ± 1.5)km/s 


HalO 


Central velocity dispersion a^ 


(8.1 ± 1.0)km/s 


HalO 


Concentration c = log(r,/rr) 


1.41 


HalO 


Core radius r^ 


16.8" 


HalO 


Half-light radius r;, 


43.8" 


HalO 



compared to the data. We summarize and discuss our results in 
Section[7l 



2. Photometry 

The photometric data for NGC 5286 was obtained from the 
Hubble Space Telescope (HST) archive. Observations were 
made with the Advanced Camera for Surveys (ACS) in the 
Wide Field Channel (WFC) with a pixel scale of 0.05 "/pixel. 
NGC 5286 is one of 65 observed clusters, being part of the 
"ACS Survey of Globular Cluster" project (Sarajedini et al.[ 



2007). The data were reduced by the "ACS Survey of Globular 



Cluster"-team as described in Sarajediru et al. (2007 1 and 
[Anderson et al.| ( [2008| l, and they also produced a star cat- 
alog, which provides a nearly complete list of all stars in 
the central 2' of the cluster, with information on position, 
V- and 1-band photometry, and some data quality parame- 
ters for each star. The Ust for NGC 5286 contains 210,729 
stars and was retrieved from the "ACS Survey of Galactic 
Globular Clusters" database, whereas the HST reference im- 
age hlsp Mcsggct Jist Mcs-wfc jigc5286 _f606w -v2 Jmg.fits in the 
V-band was downloaded from the HST- Archive. 



2.1. Center determination 

It is essential to have precise center coordinates, as the center 
influences the shape of the surface brightness profile and of the 
velocity-dispersion profile. Using the wrong center usually re- 
sults in a shallower inner surface brightness profile, and this in- 
fluences the outcome of our models. 

The center of NGC 5286 was already determined by 
Noyola & Gebhardt[ ( [2006| and [Goldsbury et~aL] ( pOTO|. With 
HST/ WFPC2 data from the F555W filter (V-band)" [Noyoir& 
IGebhardt (2006) minimize the standard deviation of star counts 



in eight segments of a circle. Goldsbury et al. (2010 1 search for 



the most symmetric point by counting stars in circular segments 
and fit ellipses to isodensity contours, using the same star catalog 
as we do. The center coordinates in right ascension a and dec- 
lination 6 are shown in Table l2] relative to our reference image. 
We apply three approaches for the center determination. 



Center determination 



a 
(h:m:s) 



Error 

(") 



:46:26.844 


-51:22:28.55 


0.5 


:46:26.831 


-51:22:27.81 


0.1 


:46:26.834 


-51:22:27.94 


0.3 


:46:26.844 


-51:22:27.88 


0.4 


:46:26.846 


-51:22:27.98 


0.4 



Noyola & Gebhardt (2006) 
Goldsbury et al.| l |2010f 
Star counts 

Cumulative distribution 
Isodensity contours 



The first method is similar to the so-called "Pie-slice 
Contours" method described by [Goldsbury et al.[ ( [2010[ ) and 



Liitzgendorf et al. ( 201 1 1. It uses star counts in opposing wedges. 



In a second approach we create a cumulative radial distribu- 
tion of the stars, as described by McLaughlin et al. (2006) and 



Liitzgendorf et al.](2011|). The crucial point of the center de- 



termination is to decide in which magnitude interval we count 
the stars, i.e. where we make the magnitude cut. This problem 
arises from the star catalog construction. In the vicinity of bright 
stars it is not possible to detect fainter stars, and the star list is 
therefore incomplete. The resulting holes in the apparent stellar 
distribution become a problem as we count the stars, and oppos- 
ing wedges seem to be asymmetric. It is therefore important to 
use only stars brighter than a certain magnitude. In addition, the 
magnitude cut must not be too low, as then there would not be 
enough stars for a statistically meaningful result and the sam- 
ple sizes would no longer be useful. The uncertainty from the 
magnitude cut is about 10 times larger than the uncertainty from 
the number of wedges. Therefore, we use a combination of the 
different parameters and weight the results with the respective 
errors. This is done for the simple star count method as well 
as for the cumulative radial distribution method. The number 
of wedges we use are between 4 and 16, and we determine the 
median and standard deviation cr„,i,ii of the resulting center co- 
ordinates. We apply this method to magnitude cuts of 19, 20, 
21, 22, and 23, and weight the resulting center coordinate with 
the respective standard deviation cr,ned to determine a center, and 
calculate the standard deviation 0-^1,, of the weighted center. As 
a result we obtain center coordinates from the star count and the 
cumulative radial distribution method and use (Jcut as a measure 
for the uncertainties. 

Alternatively, we determine the center with an approach 
that is based on the "Density Contour" method as described in 
[Goldsbury et al.| f 20 10l). The center is found by creating isoden- 
sity contours and fitting ellipses to the contours. At first we cre- 
ate a grid of 200" x 200" with a grid point every 2" , and we 
use the center determined by [Noyola & Gebhardt[ (|2006|l as our 
grid center. Assuming the center determined in Goldsbury et al. 
(2010 1 as the grid center changes the result by less than 0.04". 
At every grid point we apply a circle with a radius of 500 pixel 
(25"), find the stars inside the circle, and determine the density 
inside the circle around the grid point. Large radii of the circle 
result in smooth contour plots, but if the radius is too large, it 
reaches the edge of our star field. A contour plot is created with 
eight contour levels, and four ellipses are fit to the contours. The 
innermost and the outer three contours were ignored as they can 
be either biased by shot noise, or be incomplete. For the remain- 
ing four contours the center of the ellipse fit is determined, and 
the median is adopted as the new center, with the standard devi- 
ation of the ellipse centers as estimated uncertainty. Again, the 
choice of a magnitude cut affects the outcome of the center co- 
ordinates. Excluding faint stars smoothens the distribution, as 
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^^ 






Fig. 1. Finding chart for the different center coordinates with un 
certainties. 



Noyola & Gebhardt| ( |2006[ l: cross, [Goldsbury et al. 
( |201 0V plus sign, star count method: square, cumulative stellar 
distribution: diamond, isodensity contours: circle point. The un- 



derlying image is a synthetic I-band image (see Section 4.1 1 



fainter stars are not likely to be found in the vicinity of bright 
stars, and create under-dense regions around bright stars. This 
incompleteness is larger in the center of the cluster, as there are 
more bright stars. The difference of the center position is up to 
3" for magnitude cuts at 15 and 29. We combine the isodensity 
method for five magnitude cuts from 19 to 23. The difference 
in the center determination for these magnitude cuts is less than 
0.8", and the result therefore more robust. We calculate an error- 
weighted average, which is adopted as the new center. The stan- 
dard deviation is our estimated error, and the results are listed in 
Tabled 

Figurefllillustrates all centers from our methods with respec- 
tive errors plotted on top of a synthetic I-band ACS image (see 
Section |4.1[ i. They all lie within their error ranges, and all of 
them overlap within their errors with the uncertainty of the cen- 
ter derived by Noyola & Gebhardt (2006), and with the center 
of [Goldsbury et al. ( 2010|l. This is not u nexpected, as we use 
the same star catalog as Goldsbury et al. (|2010|l and similar ap- 
proaches to determine the center. As 'Goldsbury et al. ( 2010| l re- 
port, the center derived with the isodensity contour method is 
more reliable over a wide range of density distributions. The star 
count methods rely on the symmetry of the cluster, and depend 
more on sample size and cluster size, so we adopt the center 
calculated with the isodensity contour method. By means of the 
HST reference image we find the final position in pixel, right 
ascension a, and declination 6: 



(x„ yc) = (2934.80, 2917.59) ± (7.84, 5.88) pixel 
a = 13 : 46 : 26.847, Aa = 0.4" (J2000) 



-51 : 22 : 27.954, A(5 = 0.3" 



(1) 
(2) 



(3) 



As the distance of NGC 5286 is about 11.7kpc, the uncertainty 
of the center determination is only ~ 0.022 pc. 

2.2. Surface brightness profile 

The next step of the photometric analysis is to calculate the sur- 
face brightness profile. In the literature there are already surface 
brightness profiles of NGC 5286, one of them is from 
|Gebhardt| ( |2006l l, the other is from |Trager e"taL| ( fT995l ). 



Noyola & 



Noyola & 



14 



16 



A A A A AZ 



A^^+4%+i»j5l^ 



a 18- 



$ 20 



22- 




\ 



24 



-0.5 0.0 0.5 1.0 1.5 

log r (arcsec) 



2.0 



2.5 



Fig. 2. Comparison of the four surface brightness profiles used 
for later analysis: Our combined profile (black asterisks), the 
smoothed profile obtained from a fit with Chebychev polynomi- 
als (red plus signs), the profile from Noyola & Gebhardt (2006) l 
(cyan diamonds), and the Trager et al. '{ 1995 1 profile (green tri- 
angles, Chebychev fit on photometric points). The dotted vertical 
line denotes the core radius of 16.8" ( |Harris|1996[ l. 



Gebhardt ( 2006[ l use integrated light measurements, whereas our 
method is based on star counts in combination with integrated 
light measurement as described in Liitzgendorf et al. ( 2011| l. We 
use the HST reference image and star catalog in the V-band, and 
our center as determined in Section 2.1 via isodensity contours. 
A combination of integrated light measurement and star counts 
makes sure that the surface brightness profile is complete and 
robust. Since the star catalog does not contain every faint star, 
using only star count measurements would neglect the light con- 
tribution from faint stars. But their contribution to the light is 
taken into account in integrated light measurements. To use only 
integrated light causes problems in the determination of the sur- 
face brightness at the center. Bright stars produce shot noise and 
the resulting profile is sensitive to the choice of the bins. This 
effect is weaker with star counts, and therefore we combine both 
methods, to make our profile more robust. 

For the star counts we apply circular bins around the center, 
sum the fluxes of all stars with 14<mv<19, and divide by the 
area of each bin. The integrated light measurement is applied 
for fainter stars with mv>19, and the brighter stars are masked 
out in the reference image using a circular mask with a radius 
of 0.25" (5 HST pixels). We use the same circular bins for both 
methods and use the bi-weight to determine the distribution of 
pixel counts, and sum the fluxes of both methods. 

To obtain a photometric calibration we scale the outer part 
of our profile to the profile from |Trager et al.| ( |1995) l, using their 
Chebychev fit, and extend it with the outer points of this profile. 
In Figurel2]our adopted surface brightness profile is indicated by 
asterisks and compared to the profiles fro m [Noyola & Gebhardt] 
2006| (diamonds), and from [Trager et al. 1995 (triangles). The 



iTrager et al.[ ( [1995) ) profile is fainter and flat towards the center. 
This profile was derived from from ground based images, and 
the inner region of the cluster is not as well spatially resolved 
as with data from space. Our profile is the brightest, but due to 
the method and the data we use, the inner points of the profile 
have a high uncertainty, which comes from the Poisson statis- 
tics of the number of stars in each bin. The exact values of our 
profile are listed in Table [3] However, for later analysis we also 
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Table 3. The combined surface brightness profile in the V-band. 
AV/ and AV/, are the low and high values of the errors. 



logr 


V 


AV, 


AV„ 


[arcsec] 


mag/arcsec^ 


mag/arcsec' 


mag/arcsec^ 


-0.40 


14.83 


1.34 


0.59 


-0.19 


14.63 


0.65 


0.41 


0.00 


15.22 


0.36 


0.27 


0.18 


15.87 


0.26 


0.21 


0.48 


15.90 


0.08 


0.08 


0.65 


16.12 


0.06 


0.06 


0.88 


16.26 


0.03 


0.03 


1.00 


16.37 


0.03 


0.03 


1.10 


16.58 


0.03 


0.03 


1.18 


16.68 


0.02 


0.02 


1.30 


17.06 


0.01 


0.01 


1.40 


17.41 


0.01 


0.01 


1.48 


17.79 


0.01 


0.01 


1.60 


18.20 


0.01 


0.01 


1.70 


18.80 


0.01 


0.01 


1.78 


19.32 


0.01 


0.01 


1.85 


19.72 


0.01 


0.01 
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use a smoothed profile. This is obtained by fitting Chebychev 
polynomials to our derived profile, as indicated by plus signs in 
Figure l2] We further check whether the difference to the surface 
brightness profile from Noyola & Gebhardt ( 2006 1 comes from 
the different center position or the different methods. So we com- 
pute a profile with our method, but using the 'Noyola & Gebhardt 
( 2006| l center. At large radii the results for the different centers 
are very similar, but in the central part the difference becomes 
larger, but not more than 0.35 mag/arcsec^. The error bars are 
so large that it is difficult to say whether the choice of center in- 
fluences the shape of the profile. But the two surface brightness 
profiles with the distinct centers both tend to be higher than the 
Noyola & Gebhardt ( 2006| l profile. So the difference with their 
profile seems to come mostly from the different methods used to 
compute the surface brightness, and not from the different cen- 
ters. 



3. Spectroscopy 

3.1. Observations 

The spectroscopic data were obtained during two nights with the 
GIRAFFE spectrograph of the FLAMES instrument (Pasquini 
eral.||2002| l on 2010-05-05 and 2010-05-06 at the Very Large 
Telescope (VLT). Observations were made in ARGUS mode 
with 1:1 magnification scale, providing a total coverage of 
11. 5" X 7.3", and a sampling of 0.52" per microlens. We used the 
low resolution grating LR8, which covers the wavelength range 
from 8206 to 9400 A and has a resolving power R ^ 10,400. 
At the wavelengths 8498 A, 8542 A, and 8662 A are the Calcium 




Fig. 3. The positions of the ARGUS pointings are indicated as 
rectangles. The circle indicates the core radius of Harris ( 1996 1, 
which is 16.8". The underlying image is the synthetic I-band 
image from the HST star catalog (see Section 4.1 1. 



triplet lines, which are known to be well suited for kinematic 
analysis. Old stars, like those present in globular clusters, show 
strong features from these lines. We obtained pointings of the 
cluster's center and of adjacent regions around the core radius, 
every pointing has three exposures of 600 s each. Figure l3]illus- 
trates the arrangement of the pointings. 



3.2. Data reduction 

Data reduction includes removing cosmic rays, pipeline reduc- 
tions, and sky subtraction. We follow the procedure largely 
developed in Liitzgendorf et al. (20111, except for a change 
in the order of the steps. At the beginning, we use the 



Laplacian Cosmic Ray Identification L.A.Cosmic written by van] 
Dokkum] (2001|l to remove cosmic rays from the raw images. 
The GIRAFFE pipeline, programmed by ESO, consists of five 
recipes: In a first step, gimasterbias creates a masterbias and a 
bad pixel map from a set of raw bias frames. The recipe gimas- 
terdark corrects every raw dark frame for the bias, and creates 
a masterdark using the median. Gimasterflat creates a master- 
flat, which contains information on pixel-to-pixel variations and 
fiber-to-fiber transmissions. This recipe also locates the position 
and width of the spectra on the CCD from the flat-field lamp im- 
ages. It detects and traces the fibers on the masterflat and records 
the location and width of the detected spectra. The recipe gi- 
wavecalibration uses a ThAr arc lamp calibration frame to com- 
pute the dispersion solution and the setup specific slit geome- 
try table. The spectrum from the ThAr lamp frame is extracted 
using the fiber localization (provided by gimasterflat). The fi- 
nal pipeline recipe is giscience, which applies the calibrations 
produced by the other recipes to the cosmic ray removed object 
frames. As a result we obtain bias corrected, dark subtracted, flat 
fielded, extracted and wavelength calibrated spectra, which are 
also corrected for fiber-to-fiber transmission and pixel-to-pixel 
variations. The next step of data reduction is sky subtraction. 
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Fig. 4. (a) The synthetic image cut to the ARGUS pointings, (b) the reconstructed ARGUS image, and (c) the velocity map. Colors 
are as indicated by the color bar in units of km/s. The blue spaxels of the velocity map indicate approaching stars, the red spaxels are 
receding stars. White spaxels mean that within this spaxel no spectrum was found that could be used to determine the fit parameters. 
North is right. East is up. 



Mike Irwin developed a program, which is described in |Battaglia| 
et al.| ( [2008 1, and our program is based on this approach. It splits 
the sky spectra in continuum and sky lines, and subtracts both 
from the object spectra. We also normalize our spectra. As we 
are solely interested in the kinematic extraction from the spec- 
tra, we do not perform flux calibration. 



4. Kinematics 

The spectroscopic data contain information on the central kine- 
matics of NGC 5286, and a velocity-dispersion profile is essen- 
tial for our models. This section describes the construction of a 
velocity map, and of a velocity-dispersion profile. For the inner 
part of NGC 5286 we use the FLAMES data, for the outer part 
we use data from the Rutgers Fabry Perot (REP). 



4.1. Velocity map 

We match the reconstructed images of ARGUS produced by the 
giscience recipe to a synthetic image in the I-band, constructed 
from the photometric star catalog. Thus we identify which region 
the FLAMES data cover exactly with respect to the photometric 
center. We use the synthetic image since there is an overexposed 
region in the HST V-band image, so that the program does not 
adjust the north-eastern ARGUS images properly. At first, the 
angle between the synthetic image and each ARGUS pointing is 
evaluated, and the synthetic image is rotated by this angle. Then 
the synthetic image is blurred by convolving the image with the 
point spread function that is obtained from a Gaussian function 
with cr - 0.6", roughly corresponding to the seeing during the 
spectroscopic observations. We then cross correlate the ARGUS 
pointings with the respective blurred synthetic image. This pro- 
cedure is done iteratively, until the offset between the ARGUS 
and synthetic images is smaller than 0. 1 ARGUS pixel, which is 
achieved after four iterations at most. So, every spectrum is cor- 
related to a certain coordinate, and we obtain an area that covers 
69 X 70 spaxels. Figure |4] is an illustration of (a) the synthetic 
image, cut to match the combination of ARGUS pointings, (b) 
the reconstructed ARGUS pointings, and (c) the velocity map, 
the colors denote different values as indicated by the color bars. 
Every spectrum is a sum of individual stellar spectra, 
weighted by the luminosities of the stars. To extract the velocity 
and velocity dispersion from the integrated spectrum one needs a 
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Fig. 5. Color-magnitude diagram from the photometric HST star 
catalog of stars with low uncertainty in photometry. The star 
symbols denote the stars used for the template spectrum, dia- 
monds denote stars that contribute more than 70% of the light to 
a spaxel. 



spectral template. To find a good template we use the photomet- 
ric star catalog with information on the positions of the stars. The 
number of stars that contribute to one spaxel, and the fraction of 
light coming from a certain star in the spaxel, are calculated. 
So we can identify stars that dominate a single spaxel. Figure 
|5] shows the color-magnitude diagram from the HST catalog, 
where stars that dominate a spaxel by at least 70% are marked as 
blue diamonds. The star symbols denote seven stars, which con- 
tribute at least 84% of the light in one spaxel and are fainter than 
v(\y-\5. According to their position on the CMD, these stars are 
most probably cluster members. Those stars' spectra were com- 
bined to a master template spectrum. We found the spectral lines 
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of brighter stars are less Gaussian, possibly due to saturation, 
and therefore not suitable as template spectra. 

We use the penalized pixel-fitting (pPXF) program devel- 
oped by Cappellari & Emsellem ( 2004 1 to derive a velocity and 
velocity dispersion for each spaxel. As we have three exposures 
for every pointing and some pointings overlap, there can be sev- 
eral spectra for one spaxel. Therefore these spectra are identified, 
and sigma clipping is applied about the median of the data set for 
every spaxel, and we use the mean of the spectra for further anal- 
ysis. It is possible that only one or a few bright stars dominate the 
light in their environment. Therefore we check in which spaxel 
either a single star contributes more than 70% of the light, or less 
than five stars contribute to one spaxel. 578 spaxels of 3237 are 
affected by shot noise. Those spaxels are excluded from further 
analysis. 
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4.2. Inner velocity-dispersion profile 

To determine the velocity-dispersion profile, we follow the pro- 
cedure described in |Lutzgendorf et aL] ( |2012b} and construct cir- 
cular bins around the center of the cluster. We choose seven 
independent, non-overlapping circular bins. A higher number 
of bins results in low signal-to-noise. The spaxels that are af- 
fected by shot noise are excluded. We check different combi- 
nations of minimum star number and light contribution, and 
different angular bins, but the overall shape of the kinematic 
profile does not change much. In each circular bin the spec- 
tra are combined and we apply the pPXF method to compute 
the moments of the line-of-sight velocity distribution (LOSVD). 
As we combine all the spectra around the center in circular 
bins, we obtain the mean velocity of the circular bin {V)hin, and 

yrms - -Jy^i + ''"^' ^i'^h the rotational velocity Vrot and the ve- 
locity dispersion cr. Vrms is required for the models, and we will 
refer to the Vrms profile as velocity-dispersion profile. However, 
the Gauss-Hermite moments, which are the output of the pPXF 
program, need to be corrected as described in van der Marel & 
|Franx| ( |T993[ ). Depending on the number of fitted Gauss-Hermite 
moments, we obtain different results. We fit the first two (V, cr), 
four (V, cr, /i3, /14), or six (V, cr, /13, h^, hs, h(,) moments. 

As alternative approach we calculate the velocity-dispersion 
profile with a non-parametric fit (see Gebhardt et al. |2000b 
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Pinkney et al.||20()3l. The observed cluster spectrum is decon 



volved using the template spectrum. Deconvolution is done us- 
ing a maximum penalized likelihood (MPL) method, and one 
obtains a non-parametric LOSVD. For the fitting procedure, the 
first step is to choose an initial velocity profile in bins. This pro- 
file is convolved with the template to obtain a cluster spectrum, 
which then is used to calculate the residuals to the observed spec- 
trum. The parameters for the LOSVD are varied to obtain the 
best fit. 

The error for Vrms comes from Monte Carlo simulations, and 
is an estimate for the amount of shot noise. For every bin the 
kinematics are calculated in 1,000 realizations. The shot noise 
error is the sigma-clipped standard deviation of these 1,000 ve- 
locity dispersions. The uncertainties for the other moments of 
the LOSVD and for Vrms of the pPXF method are obtained from 
Monte Carlo simulations on the spectra. As recommended in 
Cappellari & Emsellem (2004[l, we repeat the measurement pro- 



cess for a large number of different realizations of the data. This 
is done by adding noise to the original binned cluster spectra. 
After 100 realizations we use the standard deviation of the kine- 
matic parameters as the corresponding uncertainties. 



Fig. 6. Velocity-dispersion profiles with error bars. Different 
symbols indicate different methods of calculation. M2, M4, and 
M6 are the profiles computed with pPXF, fitting only the first 
2, 4, or 6 moments, respectively. The squares denote the pro- 
file calculated with the nonparametric method. The triangles are 
from outer kinematic data, and the dotted vertical line denotes 
the core radius at 16.8" (JHams[T996|. 



Figure l6] illustrates the outcome of the different meth- 
ods. M2 is the pPXF velocity-dispersion profile when only 
the first two parameters (V, cr) are fitted, M4 when four pa- 
rameters (y, cr, /!3,/i4) are fitted, and M6 when six parameters 
(y, cr, /!3, /14, /15, /i6) are fitted. The velocity dispersion values cal- 
culated with the non-parametric method lie in between M2 and 
M4, and the shape of the velocity-dispersion profiles are all quite 
similar, but the absolute values computed with pPXF are very 
sensitive to the number of calculated moments. Therefore we 
use the non-parametric profile for further analysis. It is remark- 
able that the difference between M2, M4, and M6 is less than 1 
km/s for the innermost bin, but up to 3 km/s for the other bins. 
We tried to shift the template spectrum by the template veloc- 
ity, or use another template spectrum, but the discrepancy in the 
velocity dispersion remains. Table |4] records the results of the 
inner kinematic measurements with the non-parametric method. 
The first column denotes the radii of the bins. The following 
columns list the kinematic parameters, which are the velocities 
of each bin {V}hi„ in the cluster reference frame relative to the 
template spectrum, and the second moment Vrms with respective 
errors in units of km/s. The signal-to-noise ratio S/N is higher 
than 1 12 except for the outermost bin, where S/N=97. 

In order to estimate the radial velocity of the cluster in the 
heliocentric reference frame, we combine all spectra in the point- 
ings and measure the velocity relative to the velocity of the 
template. The velocity of the template Vrjemp results in (56.9 + 
0.6) km/s, and the velocity of all combined spectra relative to 
it is Vr,ci,ister - (2.4 + 0.6)km/s. So the heliocentric velocity of 
the cluster Vr - Vr^,emp + yr,ciusier = (59.3 ± 1.2)km/s. This is 



in agreement with the value from Harris ( 19961, Vr = (57.4 ± 
1.5) km/s. 
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Table 4. The kinematics of NGC 5286 obtained from FLAMES 
and RFP data. For the FLAMES data, {V}bin is the velocity rela- 
tive to the template velocity. 

FLAMES DATA 



r (V)b,„ A(V),„„ 

[arcsec] [km/s] [km/s] 



0.52 
1.56 
3.12 
5.46 
8.06 
12.48 
20.54 



[arcsec] 

37.50 
60.00 



V„„, AV„„, 
[km/s] [km/s] 



2.4 
4.6 

2.7 
2.1 
2.7 
2.4 
2.1 



0.4 
0.4 
0.3 
0.4 
0.3 
0.4 
0.4 



9.1 

10.4 
9.5 
10.1 
9.3 
8.8 
8.2 



2.2 
1.6 
1.0 
0.8 
0.7 
0.7 
0.7 



RFP DATA 



[km/s] 

59.0 
57.9 



[km/s] 

0.8 
0.8 



' rnis 

[km/s] 

7.7 
7.3 



AV„„, 
[km/s] 

0.5 
0.6 



S/N 



112 
133 
140 

142 
138 
134 
97 



Number 
of stars 

109 

80 




Fig. 7. The mean radial velocities with error bars of stars at a 
distance r > 25" plotted versus position angle (p for 6 segments 
from RFP data. The solid line is the fit to the data. 

velocity-dispersion profile. We obtain cr^ = (9.3 + 0.4) km/s, this 
is higher than the central velocity dispersion of (8.1 ± 1.0) km/s 
from Harris ( 1996 1, but the values agree within their uncertainty 
limits. 



4.3. Outer velocity-dispersion profile 

In addition to spectroscopic data in the center we use a second 
data set for larger radii, i.e. up to r~75". The data were ob- 
tained with the 4-m Victor M. Blanco Telescope at the Cerro 
Tololo Inter-American Observatory (CTIO). Observations were 
made from 1994-05-31 to 1994-06-03, with the Rutgers Fabry 
Perot (RFP). The used methods of observation and data reduc- 
tion are similar to those described in |Gebhardt e t al.|(|1997[l for 
M15. The data set for NGC 5286 contains 1,165 velocities of 
individual stars with errors, flux, and information about their po- 
sition relative to a center in arcsec. With an RFP image of the 
pointing, the HST image, and the iraf recipes daofind, ccmap, 
and xyxymatch, we transform the coordinates to the world coor- 
dinate system (WCS). 

To determine the velocity dispersion we draw circular bins 
around the center and use only stars with low error in veloc- 
ity and with a high flux. This is done to make sure that only 
stars with a well-determined velocity are used for further analy- 
sis. Crowding in the cluster's center is known to cause problems 
in the determination of velocity. Background light contaminates 
the measurements for fainter stars, and this biases the velocities 
of the stars to the mean velocity of the cluster. Therefore, the ve- 
locities obtained from Fabry-Perot measurements are often too 
low. We use only stars in the outer part of the cluster, where the 
crowding efifect is less relevant. With the maximum likelihood 
method introduced by |Pryor & Meylan| (|1993|) we calculate the 
velocity dispersion Vnns and mean velocity (V)/,,,,. The obtained 
values are listed in Table HI The last column displays the num- 
ber of stars that were used to determine Vrms for each bin. Both 
values of {V}bi„ are in agreement with the heliocentric velocity 
of the cluster from Harris ( 1996 ), Vr - (57.4 ± 1.5) km/s. 

We also compute the efifective velocity dispersion cTe, as de- 
fined in |Gultekin etaL] ( [2009l ) 



cr? = 



//' Vlj(r)dr 



J^/I(r)dr 



(4) 



where R,, is the effective radius or half-light radius (43.8'' 
|19961 l, I(r) is the surface brightness profile, and V,.„, 



Harris 



4.4. Rotation 

Some globular clusters are known to rotate, for example 
wCentauri and 47 Tucanae ( [Meylan & Mayor||1986| ). As NGC 
5286 is not totally spherical and has an ellipticity e - 0.12, com- 
parable to the ellipticity of w Centauri (e = 0. 17) and 47 Tucanae 
(e = 0.09, aU eUipticity values are from [ White & Shawl|1987) , 
we investigate a possible rotation of NGC 5286 with our two 
data sets. For the outer part of the cluster we use the RFP data. 
We separate the cluster in 6 to 18 segments and find all stars 
with a distance from the cluster center of at least 25". The mean 
velocity of the stars in these areas is determined with the maxi- 
mum likelihood method. In Figure [Tithe mean radial velocity V 
of the stars in the segments is plotted against the position angle 
(p. For six segments there are more than 85 stars in each bin, for 
18 segments at least 20 stars. We fit the function / 



/ = A sin((/> + t//) + Vr, 



(5) 



the 



with the parameters amplitude A, phase angle i/f, and heliocen- 
tric cluster velocity Vr to the data points for all segment numbers 
from 6 to 18. The best-fit parameters are averaged, and we ob- 
tain for the mean heliocentric velocity Vr = (58.3 + 1.1) km/s, 
and for the amplitude A = (2.3 ± 1.5) km/s. The phase t// is 
(-0.01 ±0.7)rad, and the position angle of the rotation axis is 
(89 ± 40)°, measured positive from north to east. White & Shawl 
( 1987[ l found the orientation of the major axis at 0=21°. 

The spectroscopic FLAMES data is used to investigate in- 
dications of rotation in the inner part of the cluster. In some 
segments the spectra have low signal-to-noise and therefore the 
mean radial velocities have large error bars of up to more than 
3 km/s. For the sine fit we use only segments with a signal-to- 
noise higher than 40. However, the result for the phase angle if/ 
depends highly on the number of fitted segments, and therefore 
we find no significant indication for rotation in the inner 25.5"of 
NGC 5286. 

5. Jeans models 

Jeans models provide an analytical approach to fit kinematics of 
the globular cluster NGC 5286, and we use an IDL routine devel- 
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Fig. 8. Different spherical Jeans models fitted to our combined 
surface brightness profile for fixed anisotropy j6=0. 1. The up- 
per panel shows the models together with the data for differ- 
ent black-hole masses, the thick black line indicates the best fit. 
The lower panel shows the x^ values as a function of black-hole 
mass, and the dashed line denotes Ay^ - 1 . The best fit with x^ 
- 1.89 requires a black-hole mass of M. = (1.5+)°) ■ 10^ Mq. 



oped bylCappellari ( |2008 1. In order to solve the Jeans equations, 
the mass density q is required as an input. We use the surface 



brightness profile for NGC 5286 derived in Section 2.2 and de 



project it to obtain the luminosity density j of the cluster, which 
is related to the mass density q by the mass-to-light ratio (M/L). 
The surface brightness profile is deprojected using the Multi- 
Gaussian expansion (MGE) method, developed by |Emsellem| 
etal.|([T994), using the IDL routine written by|Cappellari (2002). 



With the MGE formalism there are Jeans equations for every sin- 
gle Gaussian component, and after solving the equations, we ob- 
tain the second velocity moment f^„(R) along the line-of-sight. 
The results are fitted to the observed velocity-dispersion data 



points, as computed in Sections 4.2 and 4.3 It is possible to vary 
the values of anisotropy fi-\ - cr'-gla-j^ and M/L for every in- 
dividual MGE Gaussian. The velocity-dispersion profile can be 
scaled to the observations by adjusting M/L. Therefore all mod- 
els meet in one point at which the model is scaled. 



5. 7. Spherical Jeans models with constant M/L 

The simplest model is a spherical Jeans model. We consider 
different values of anisotropy /3, and different surface bright- 
ness profiles. The profiles used are the Chebychev fit from 
Trager et al."1[T995] [Noyola & Gebhardt|[2006l our derived pro- 



file from a combination of star counts and integrated light mea 
surements (see Table[3j combined profile), and the profile from a 
fit of Chebychev polynomials to the combined profile to make it 
smooth (smooth profile). The best fits are obtained with our com- 
bined profile. For each profile we test different values of constant 
anisotropy /? in the interval of [-0.2, 0.3] over the entire cluster 
radius. We fit a global M/L, and the best fit is obtained with ra- 
dial anisotropy /3-0.1, M.=L5-10^ Mq, and constant M/L=L07 
(all values of M/L are in solar units of Mq/Lvq). The outcome 
is shown in Figure^ For the other profiles the best fit values 
are listed in Table ^with errors from the Ax^-l limit. iTrager 
et al. ( I995|l do not provide uncertainty limits to their surface 
brightness profile, therefore the uncertainty in M,o, comes only 



from the uncertainty in M/L. The value from Harris ( 1996 1 for 
the integrated V-band magnitude is 7.34, this corresponds to a 
total luminosity L,„, of 2.40 ■ 10^ Lq. The [Trager et al!] ( |l995 i 
profile has too low resolution at the center and is therefore not 
as suitable for our analysis as the other profiles. 

We also assume a radial shape of anisotropy /? oc r, with 
isotropy in the very center of the cluster, and rising anisotropy 
towards the outskirts. At the outermost Gaussian the anisotropy 
j3oui is highest. This approach is justified as the relaxation tim e is 
short in the cluster's center (treia.x ~ 2.5 10^ yr, Harris 1996 1, and 
higher at larger radii (at the half- mass radius ti-eiax ~ 1-3 ■ 10^ yr, 
Harris 1996) . |Liitzgendorf et al.| ( |20I 1) ran an A^-body simula- 
tion of a globular cluster with initial anisotropy and showed that 
it becomes more isotropic after few relaxation time s. Assuming 
an age of 11.7 ■10'' yr (Rak os & Schombert|[2005l ), NGC 5286 
is already 46 relaxation times old in its center, and about 9 
relaxation times at the half-mass radius. Therefore the cluster 
should be isotropic in the central part at least. In this case, the 
fits have similar values of x^ and obtain similar results for all 
values of Pom at a given surface brightness profile. Our com- 
bined profile finds values for the black-hole mass in the range 
of (1.9-2.3)-10^Mq and M/L=1.07-1.09 at all values of /?. 
Figure^ is a contour plot of ;t'^ as a function of j6 and M, from 
the combined surface brightness profile. For these models we 
use a constant and global M/L = 1.1, and do not fit the best 
value of M/L, since we are interested in showing the influence 
of yS and M. on the quality of the fit. If M/L is fitted, it usually 
has a value around 1 . 1 anyway. In the case of radially varying 
anisotropy the abscissa denotes the value of y8o,„, the anisotropy 
in the outermost part of the cluster. The value of x^ is labeled 
on the contours. For constant /J the area with the best models is 
smaller than foryS oc r, where fits over a large range of /3 are of 
similar goodness. We conclude that anisotropy in the outer part 
of the cluster has no big effect on the resulting black-hole mass 
in Jeans models as expected, but the central value of anisotropy 
is crucial. 

From the Ax^ = 1 limit we obtain an estimate of the un- 
certainty of a single model. However, to estimate the error of 
M/L and M, we also run Monte Carlo simulations. In 1,000 
runs we vary the shape of the combined surface brightness pro- 
file to estimate the effect on the outcome of the Jeans models. 
We change only the six innermost points of the profile, as they 
have the largest uncertainties and are crucial to determine the 
mass of a central black hole. This is done for /? in the range 
of -0.20 to 0.20. From the 68% confidence limit we obtain 
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Table 5. The spherical Jeans model best-fit values with constant M/L. 

surface brightness profile /? M/L M, M,„, L,g, x^ 

constant [Mq/Lv.o] [lO^ Mq] [iO'^Mq] [iO^Lq] 



combined profile 
smooth profile 



Noyola & Gebhardt 1 2006 1 



Trager et al. ( 1995 i 



0.10 
0.05 
0.10 
0.15 



1.07 ±0.03 



1 07+°-"'* 
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Fig. 9. ;i'^ contours describing the agreement between the data and spherical Jeans models of the combined profile with constant 
mass-to-light ratio M/L=l.l, the numbers denote the respective value of ;(f^. Left: /3 - const., right: p oc r. In the case of radially 
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Fig. 10. The result of Monte Carlo simulations on the surface brightness profile. After 1000 realizations for every value of constant 
y6 the dots denote the median values, the contours are the 68% (dotted line) and 95% (dash-dotted line) confidence limits. The left 
panel displays the mass-to-light ratio M/Ly, the right panel is the black-hole mass. 



the uncertainties SMfLy = 0.01 and 6M. = 0.4- IO^Mq. For 
a confidence limit of 95% we find (JM/Ly = 0.08 and SM, - 
1.2 ■ lO"* Mq. By varying both, the surface brightness profile and 
the velocity-dispersion profile, we run another Monte Carlo sim- 
ulation and find higher uncertainties of6M.fLy - 0.04 and 5M, - 
1 .0 ■ 10^ Mo for the 68% confidence limit, and cJM/Ly = 0. 16 and 



SM. - 4.0 ■ 10^ Mq with 95% confidence. The difference of the 
uncertainties by varying both profiles or only the surface bright- 
ness profile give information on the contribution from the dif- 
ferent profiles. The uncertainty of M/L seems to come mostly 
from the uncertainty of the kinematic profile, and the influence 
on the black-hole mass is also dominated by the kinematic pro- 
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file. Nevertheless, one must not underestimate the importance of 
the shape of the surface brightness profile. The Ay^ - 1 limit of 
a single model gives also uncertainties in the same order of mag- 
nitude as the 68% confidence limits, though some models have a 
larger range of uncertainty than from Monte Carlo simulations. 



Figure 10 shows the median values of these 1,000 realizations of 



the surface brightness profile as red points and as a function of 
y6, and the 68% and 95% confidence limits are the light grey and 
grey shaded contours. The uncertainty of M/L is plotted in the 
left panel, M, is plotted in the right panel. 

All our results come from Jeans models in which we neglect 
the PSF effects as described in Appendix A of Cappellari (2008]l. 
To test the importance of this effects we run spherical Jeans mod- 
els and convolve the model velocity dispersion with the instru- 
mental PSF before fitting to the observed data points. We use 
the combined surface brightness profile and constant values of 
anisotropy. With a seeing of 0.6"and a pixel scale of 0.52 "/pixel 
the values for black-hole mass, total cluster mass, and mass-to 
light ratio decrease by less than 0.4%. As the convolution is time 
consuming but has no big effect on the result, we disregard it in 
all further Jeans models. 



5.2. Spherical Jeans models with a varying M/L profile 

A'-body simulations of globular clusters usually result in mass- 
to-light ratios that vary with radius. This is expected from mass 
segregation of the stars, as more massive stars, in particular dark 
remnants, wander to the center of the cluster, and low-mass stars 
to the outskirts. We want to investigate the effect of a varying 
M/L profile on the Jeans models, and therefore calculate it from 
A^-body simulations, and use it as input to the Jeans models. The 

it comes from the best fit of 



profile is displayed in Figure 11 



a grid of A^-body simulations, not including the effect of a tidal 
field (see Section 6.1 1, to our profile from a combination of star 
counts and integrated light. This profile is scaled to the velocity- 
dispersion profile in the Jeans models to obtain absolute values 
for M/L in the V-band. We use only the combined profile, as it 
gives the best fit for constant M/L, and this profile was also used 
to determine the shape of the M/L profile. 

With a constant value of /3 - 0, the best-fit black-hole mass 
is (2.0 +[|^) ■ 10^ Mo with;\^2 = 2.21 and the mean value of M/L 
is ~ 1.5. A better fit is obtained with a constant anisotropy of y6 
= 0.15 with (0.9 t_oi) ■ 10^ Mo, M,o, = (4.1 + 0.2) ■ 10^ Mo and;^'^ 
= 0.99. Except for the total cluster mass, the results are in agree- 
ment with the spherical Jeans models with constant M/L, they 
yield similar masses for a central black hole. The higher clus- 
ter mass comes from the M/L profile, which has higher values 
of M/L at small radii, where our surface brightness profile has 
the largest error bars. It also rises towards the outskirts, but our 
velocity-dispersion measurements extend only to 75". Therefore 
the variable M/L profile adds mass to the cluster without having 
a big effect on the outcome of the Jeans models for M,. 

5.3. Axisymmetric Jeans models with constant M/L 

The axisymmetric model requires the 2D surface brightness pro- 
file, which we obtain with the MGE-fitsectors IDL package 
written by Cappellari (20021. We use the HST image and three 



2MASS images in the J-band, which were combined to one im- 
age using the tool Montage. The routine find-galaxy.pm mea- 
sures the ellipticity e - \ - b/a from a two dimensional image. 
One parameter of find-galaxy.pro is the fraction of image pixels 
that are used to measure the outcome. We use different fractions 




10 
Radius [arcsec] 



Fig. 11. Mass-to-light ratio as a function of radius, calculated 
from the best fitting A^-body model. 



(0.1, 0.15, 0.2, 0.25, and 0.3), and compute the mean values for 
ellipticity, which is 0.12 + 0.02 in the HST, and 0.14 ±0.04 in 
the 2MASS image. 



Both values are in agreement with the value from White & 



Shawl ( 1987) of e - 0.12. The routine sectors photometry. pro 
determines a 2D surface brightness profile from the HST and 
2MASS images. The 2MASS photometry is matched to the 
HST flux in the range from 2" to 100", and scaled to the 
HST flux. The HST measurements are further used for the 
central part of the cluster up to 60" , the 2MASS photometry 
for the outer part (13"- 400"). The photometry is then scaled 
to the [Trager eTal] ( [T995] l profile at 25"- 158". The routine 
MGEJit-sectors.pro fits a 2D MGE model to the measurements, 
and to every Gaussian of the fit there is an axial ratio Q < qi^ <\ 
assigned. The total luminosity is L,o,=2.67 ■ 10^ Lo. Varying the 
range where the photometry is scaled to |Trager et al.| ( [T995| l in- 
troduces an uncertainty of 9 ■ lO"* Lo. 

Axisymmetric Jeans models have an additional input pa- 
rameter, the inclination angle /, which is the angle between 
the rotation axis and the line-of-sight. An inclination angle 
/ = 90° means the cluster is seen edge-on, and / = 0° corresponds 
to face-on. But in the face-on case we would probably see no 
indication of rotation, as we do in the outer part of the cluster, 
and the ellipticity e would be zero. We run axisymmetric Jeans 
models for / = 45°, 67°, and 90°, and constant/? in the interval 
of [-0.2, 0.3]. M/L is assumed to be constant and fitted in every 
model. 

The best fits at every value of / are listed in Table l6] All 
errors are from the Ax^-\ limit, the uncertainty in M,o, comes 
from the uncertainty in M/L and L,„,. At lower inclination we 
obtain better fits, and in contrast to spherical Jeans models, the 
best fit is at tangential anisotropy (fi - -0.1). At higher values of 
i the isotropic case fits best. But all models obtain similar results 
for M/L, M., and M,o,, which agree within their uncertainties, 
and are also consistent with the results from the spherical Jeans 
models. 



6. A^-body simulations 

6.1. Fitting the data to a grid of N -body simulations 

In addition to Jeans models we use a grid of A^-body models 
and fit the models to the data of NGC 5286. The simulations 
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Table 6. The axisymmetric Jeans model best-fit values. 



i p MIL M. M,„, 

n constant [Mq/Lv.g] [iO^Mg] [iO^Mq] 



^tot 

10= Lo] 



45 


-0.10 


i.08!°:«^ 


1 5+'* 


67 


0.00 


1.04!!;- 


1 9+"* 


90 


0.00 


1.04!;;- 


1 9+''* 



X^ 



1 + 0.13 2.67 ±0.09 1.16 



2.76 ±0.13 2.67 ±0.09 1.36 



2.76 ±0.13 2.67 ±0.09 1.39 



are described in [McNamara et al.| ( |2012| l. All simulations use 



the code NBO DY6 (|Aarsethl|2003| and start from King ((T962jl 
models with a Kroupa (2001) initial mass function between 0. 1 
and 100 Mq. The stars evolve according to the stellar evolution 
model of [Hurley et al.| ( [2000| l. The simulations run for 12 Gyr 
and we use 10 snapshots with a 50 Myr interval in between, start- 
ing at 1 1 Gyr. We overlay the snapshots for comparison with the 
data. The cluster is then scaled up to NGC 5286 using the scal- 
ing as described by Jalali et al. ( 2012[ ). To obtain smooth surface 
brightness and velocity-dispersion profiles from the simulations 
we use the infinite projection method of jMashchenko & Sills| 
(J2005 ), which averages the stellar positions over all possible ro- 
tations of the cluster relative to the observer. 

Simulations cover a three-dimensional grid with different 
values of black-hole mass, concentration c - log(r,/ff), and ini- 
tial half mass radius R///, and we search within this grid to obtain 
the model which best fits the observational data. The simula- 
tions are run without an external tidal field, so cluster evolution 
is driven by two-body relaxation and stellar evolution. For the 
fit to the A^-body model we use our combined surface brightness 
profile, as it also results in the lowest x^ with spherical Jeans 
models. The error bars in the outskirt of the profile are less than 
1%, which gives too much weight on the fit of the surface den- 



sity. Therefore we use the approach of McLaughlin & van der 
Marel|([2005) to estimate the error ctj of the Trager et al. '{ 1995|l 



profile from the relative weights of each data point (w,) and a 
constant tr^, with cr,- ^cr^/w,. The weights were estimated by 
Trager et aLl ( |1995| l, and [McLaughUn & van der M arel (200^ 
give the value of cr^ =0.128. We use these errors for the com- 
bined surface brightness profile, starting at ~ 3". 

The best fitting models for no IMBH, an IMBH with 1%, 



and with 2% of the total cluster mass are displayed in Figure 12 



The surface density is shown in the left panel, and the velocity 
dispersion in the right panel. We convolved the model profiles to 
take the effects of seeing and pixel size into account. The change 
of the model profile is small except for the center where R< 0.5", 
beyond our data points. All models have similar results for the 
surface brightness at radii R > 5" and differ only at the center. 
Regarding only the surface brightness profile, the 2% IMBH 
model fits better than the 1% IMBH or the no-IMBH model. 
However, the velocity-dispersion profile of the 2% IMBH model 
rises too much in the center of the cluster, and does not fit the 
innermost data point. The 1 % model also rises towards the cen- 
ter, but not as much. The best overall fit contains an IMBH with 
M./M,o, = 0.9% and;^^ = 11.3 (xl^^^Q39). This model gives a 
total cluster mass of M,o, = (4.38 + 0. 18) ■ 10^ Mq, therefore the 
IMBH has M. = (3.9 ± 2.0) ■ 10^ Mq. The best-fitting no-IMBH 
model has M^, =4.48 ■ 10^ Mq with;t'2 = 13.4. This case can be 
excluded at the 68% confidence limit. 

Figure 13 shows the;i'^ contour lines as a function of the start 
parameters of the simulations for a 1 % IMBH model. The lowest 



X^ for the surface density are at a region of models which start at 
a concentration of about c = 2 (panel a). The value of the initial 
half-mass radius R/„ is not much constrained. For the velocity 
dispersion (panel b), best fits are in a narrow band going from 
Rhi - 4.5 pc and low concentration to higher concentration mod- 
els and R//, = 2 pc. McNamara et al. ( 2012| l observe the same for 
NGC 6266 and argue that this is due to the dynamical evolution 
and the performed scaling of the cluster. However, the combined 
reduced x^^ considering the fit to both profiles (panel c) comes 
from a region around c = 1 .6 and R/f, =2.1 pc and has the value 

;rL=o.39. 



6.2. Testing the spherical Jeans models with an N-body 
simulation 

In order to test whether Jeans models give reliable results for the 
total cluster mass, we use the A^-body simulation which is clos- 
est to the best-fitting one and analyze it in more detail in this 
section. This simulation started with an IMBH mass of 1%, a 
half-mass radius is R/// = 2.0 pc, and an initial concentration of 
c = 1.5. The simulation leads to a cluster mass of 4.94 ■ 10^ Mg. 
We compute the surface brightness and velocity-dispersion pro- 
files for this simulation and perform a spherical Jeans analy- 
sis with these profiles. To compute the velocity dispersion we 
make three perpendicular projections and compute a velocity- 
dispersion profile for each projection with the method of |Pryor| 



|& Meylan ( 1993 i. The three projections are averaged to a final 
velocity-dispersion profile. We also vary the limiting magnitude 
of stars that contribute to the velocity-dispersion profile. Fainter 
stars tend to be at larger radii due to mass segregation, and there- 
fore increase the velocity-dispersion profile there. Computing 
the profile with brighter stars only underestimates the velocity 
dispersion at large radii, and the profile, which rises to the cen- 
ter, is too steep. Radial anisotropy also steepens the velocity- 
dispersion profile. 

We run Jeans models with constant M/L and the radially 
varying M/L profile from the simulation, and different constant 
values of j3. The limiting magnitude was chosen to be either 
mv= 19, 20, 22, 24, 26, 28, or 30. The influence of varying the 
limiting magnitudes, the anisotropy parameter /3, and the M/L 
profile on the total cluster mass is shown in Figure 14 The ex- 



tracted velocity dispersion rises as a function of limiting mag- 
nitude. This results in a rise of total cluster mass as a function 
of limiting magnitude. We run these tests with Jeans models of 
constant M/L as well as for radially varying M/L and find the 
same behavior for both. However, the recovered M,„, is system- 
atically lower for constant M/L, as constant M/L underestimates 
the cluster mass at large radii. Also the value of /3 influences the 
total cluster mass M,o, found by the Jeans models. Tangential 
anisotropy results in higher M,,,, than radial anisotropy. The ra- 
dially varying M/L profile at a magnitude cut of mi/=26 andy6=0 
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Fig, 12. Best fitting A^-body model without tidal field and with outer kinematic data. The surface density is plotted on the left panel, 
the velocity dispersion on the right panel. The different colors denote the best models for a black hole with 1 % of the cluster mass 
(red), 2% of the cluster mass (green), and for no black hole (blue). The data points are denoted as black symbols. 
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Fig. 13. ;f^ contour lines as function of start parameters King concentration c and initial half-mass radius Rhi for the surface density 
(a), velocity dispersion (b), and both profiles combined (c). The red star indicates the best fit. 



obtains M,o, - 4.96 ■ 10^ Mq, which is almost perfectly matching 
the input value of 4.94 ■ 10^ Mq. As our spectral observations do 
not contain such faint stars, our measured velocity dispersion for 
NGC 5286 is probably also too low. This can be the reason why 
the Jeans models obtain a lower total cluster mass than the A^- 
body simulation, even after taking the radially varying M/L into 
account. 



7. Conclusions 

7.1. Summary 

Using photometric and kinematic data, we modeled the glob- 
ular cluster NGC 5286 and found indications for the presence 
of a central black hole. Our photometric data set consists of an 
ACS/HST image, and the star catalog constructed by 'A nderson] 
let al. (2008) from ACS/HST observations. We used these to de- 
termine the photometric center with three different methods and 
confirmed the results of |Goldsbury et al.] ( |2010| l. Further, we pro- 



duced a surface brightness profile from a combination of star 
counts and integrated light measurements. We used two kine- 
matic data sets, for the inner part of the cluster up to ~ 25" we 
have ground based spectroscopic data from the integral field unit 
ARGUS at the VLT This data set was reduced and used to pro- 
duce a velocity map, and a velocity-dispersion profile. The sec- 
ond data set was obtained with the Rutgers Fabry Perot at the 
CTIO and contains velocity measurements of single stars. We 
used the outer part up to a radius of ~ 75" to complement our 
velocity-dispersion profile. Examining both data sets for indi- 
cation of rotation, we found rotation only in the region beyond 
25" with a rotation velocity of (2.3 ± 1.5)km/s and a rotation 
axis consistent with the minor axis of the cluster. 

With the information provided by the surface brightness pro- 
file we ran Jeans models and compared the results for the veloc- 
ity dispersion to our profile, x^ statistics were used to find the 
best fit. For the spherical models we used our surface bright- 
ness profile as well as the profiles from Trager et al. ( 1995| l, 
and [Noyola & Gebhardt| ( |2006| l. With different central black- 
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Fig. 14. Total cluster mass M,o, in Mq from spherical Jeans models, performed on the profiles of an A^-body simulation. For the 
surface brightness profile we use different limiting magnitudes my, thus changing the number of stars contributing to the profile. 
Left panel: constant M/L. Right panel: radially varying M/L profile. The solid black line is the input value of M,,,, in the simulation. 



hole mass and anisotropy behavior the outcome of the model 
was scaled to our velocity-dispersion profile to determine the 
value of a constant mass-to-light ratio. Further, we used an M/L 
profile computed from A^-body simulations as input. As a sec- 
ond approach we computed a two-dimensional surface bright- 
ness profile from the ACS/HST image and 2MASS images and 
ran axisymmetric Jeans models with different black-hole masses, 
anisotropy, and inclination. We also fit the surface brightness 
profile and the velocity-dispersion profile from a grid of A^- 
body models without tidal field to our data. All our models re- 
quire a central black hole for their best fit. Jeans models find a 
best fitting black-hole mass of about (1.5±1.0)- IO-'Mq, while 
N-body modes require a black hole of 0.9% of the total clus- 
ter mass, corresponding to M, = (3.9+2.0) ■ 10^ Mq. The errors 
are the 68% confidence limits. The models also provide the to- 
tal mass of the cluster, and the ratio of the black-hole mass 
to the total cluster mass is lower for Jeans models, M.fM,,,, ~ 
0.5% - 0.7%. A^-body simulations obtain a total cluster mass 
M,ot - (4.38 ± 0.18) ■ 10^ Mq. The total cluster mass of spheri- 
cal Jeans models is up to 34% lower, depending on the value 
of M/L. For constant M/L we obtain a value of 1 .07 M/Ly and 
Mfo, = (2.89 + 0.10)- IO^Mq, but with the varying M/L pro- 
file from A^-body simulations the cluster mass is higher (M,o,= 
(4.1±O.2)-1O^''M0). 



7.2. Discussion 

The value of M/L from our Jeans models with constant M/L 
is about I.O7M0/LV0. Both spherical and axisymmetric Jeans 
models find similar values. The surface brightness profiles from 
Noyola & Gebhardt] (J2006 ) and [Trager et aL] ( |1995[ ) result in 
higher values of 1.11 Mq/Lv_q and 1 .22 Mq/Lv_q, respectively. 
This is probably due to the lower value of total luminosity with 
these profiles (see Table B). The integrated V-band magnitude 
of |Harris| ( |1996[ ) is 7.34, corresponding to a total luminosity of 



2.7 ■ 10^ Lq. This value is in good agreement with the total lumi- 
nosity from our profiles. 



Pryor & Meylan ( 1993 1 found a value of M/L = 2. 1 Mq/L 



v.o 



from integrated light measurements and isotropic King mod- 



els. McLaughlin & van der Marel 
et al. ( 1995|l profile and data from ' 



( 2005 1 used the Trager 



iarris (19961 to compute 



M/L ratios. They used a population-synthesis model (M/L = 
(1.87 ±0.16) Mq/Li/q) and dynamical models to fit M/L (King 
fit: 0.99 +;}49^ Wilson fit: 0.94 +°-;^^, and Power-law fit: 0.91 
-035 ™ solar units, respectively), and obtained overall lower 
values than Pryor & Meylan ( |1993[ l, in better agreement with 
our Jeans model with constant M/L of ~1.07 Mq/Lvo . 'Marastonl 
(2005) also calculated M/L ratios from evolutionary population 
synthesis models. The results of the Sloan Digital Sky Survey 
(SDSS) r-band filter at 622 nm can be compared to our values 
of M/Ly in the V-band at 606 nm. With the |Kroupa] ( |2001| l ini- 
tial mass function, a metallicity of [Z/H] = -1.35 and an age of 
11-12 Gyr, their model gives higher values of M/L= 1.87-2. 02. 
Our M/L profile from A^-body simulations has a global value 
of M/Lv of about (1.75 ±0.06) Mo/L^q after scaling to the 
velocity-dispersion profile, which is in rough agreement with the 
M/L ratio of the population synthesis models. 

The differences in the total cluster mass between our mod- 
els can be explained by the influence of the M/L profile. 
With the M/L profile from A^-body simulations as input to 



our spherical Jeans models in Section 5.2 we obtain a higher 
cluster mass of M,o, = (4. 1 ± 0.2) • 10^ Mq, in better agreement 
with the value found from the A^-body simulations of M,^, = 
(4.38 ± 0.18) ■ 10^ Mq. Also the Jeans models that we ran on the 
profiles of an A^-body simulations in Section |6.2| confirm the 
influence of a radially varying M/L profile on the total cluster 
mass. In the range of 2"< r < 80" the M/L profile is lower than 
average. The extra mass, compared to a constant M/L, is mostly 
at larger radii. Since we have no velocity measurements for r 
>75", the additional mass at large radii does not significantly 
influence our velocity-dispersion measurements. A rise in M/L 
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towards the outer cluster parts could also explain why the dy- 
namical M/L values of McLaughlin & van der Marel ( 2005 | l 
are smaller than the population synthesis ones. This stresses 
the importance of obtaining velocity-dispersion measurements 
in the outer cluster parts when trying to determine the total clus- 
ter mass. Furthermore, the velocity-dispersion measurements are 
based on bright stars. Those stars are more massive than average, 
but have lower velocities. This underestimates the velocity dis- 



persion and the total cluster mass, as we showed in Section 6.2 



For the determination of a central black-hole mass, the value 
of the anisotropy fi is important. Anisotropy in the outer part 
of the cluster does not change the result much, but a radial 
anisotropy in the center can easily be misinterpreted as a cen- 
tral black hole, whereas unknown tangential anisotropy results 
in an underestimation of the black-hole mass. Assuming ra- 
dial anisotropy of /3- 0.30 in the center, a central black hole 
is not necessary to explain the velocity-dispersion profile of 
NGC 5286. To measure the amount of anisotropy, proper motion 
measurements are needed in addition to radial velocity disper- 
sion. However, proper motion measurements require a long time 
baseline and high accuracy. Crowding in the center causes fur- 
ther problems. Therefore the only reliable method to constrain 
the anisotropy are A^-body simulations. Liitzgendorf et al. ( 2011| l 
ran a simulation of a highly anisotropic cluster and showed that 
the central part of a cluster becomes isotropic after a few relax- 
ation times. As NGC 5286 is already 9 relaxation times old at 
its half-mass radius, we expect not much anisotropy (^| < 0.2). 
Especially in its central part, where the relaxation is even faster, 
the cluster should be isotropic. Nevertheless our Jeans models 
obtain better fits for radial anisotropy than for isotropy. Mass 
segregation could be one reason to explain this: The low-mass 
stars, which are mostly in the outer part of the cluster, are too 
faint for detection, and the velocity-dispersion profile at larger 
radii can therefore be underestimated compared to inner radii. 
Our velocity-dispersion profile may be too steep, and the Jeans 
models misinterpret this as radial anisotropy. 

The Jeans models find smaller values for the mass of the 
black hole than the A^-body models, but both indicate a black 
hole with less than 1 % of the total cluster mass. This upper limit 
is higher than the black hole mass fraction usually measured for 
galaxies. Haring & Rix ( 2004) found that the fraction of a central 
massive dark object in nearby galaxies is around 0.14% + 0.04% 
of the bulge mass of the host system. 

We have 1-cr detections for an IMBH, so the no-black-hole 
case cannot be excluded. A central black hole should not be 
more massive than 6.0 ■ 10"^ Mq, as for A^-body simulations this 
is where A^^ - 5. Jeans models find an upper limit of M, = 
5.5 ■ lO"* Mq with 95% confidence level from Monte Carlo simu- 
lations. But independent of the method we use, all our methods 
prefer the case with an IMBH over the no-black-hole case. The 
result from A^ -body simulations is M,=(3.9 ± 2.0) ■ 10^ Mq, from 
Jeans models the black-hole mass is about 50% lower. 

Our modeling shows that the derived uncertainties for the 
kinematic profile are too high to precisely constrain the black- 
hole mass and anisotropy. In the innermost bin the uncertainty 
is more than 2.2 km/s, which is 24% of the measured value. 
However, these high error bars are realistic, as the determina- 
tion of the velocity dispersion currently shows some ambigu- 
ity. We used two different methods and varied the parameters 
of the pPXF method. Thus we found scatter of up to 3 km/s. 
Further investigation of this inconsistency is therefore needed 
to obtain better constrained results. Alternatively, observations 
with higher spectral resolution or of individual stars could also 
clear up this issue. 



NGC 5286 is an interesting case for dynamical modeling 
and could host an intermediate-mass black hole. Our research 
group has IFU data for six more globular clusters in the south- 
ern hemisphere (Liitzgendorf et al. 2012a). These data are used 
to constrain possible IMBHs inside these clusters. This infor- 
mation can be used to decide which other clusters should be 
examined, also in the northern hemisphere. Located between 
galaxies and globular clusters, at cr,, ~ 15 - 70 km/s, are systems 
like dwarf galaxies, which could also contain intermediate-mass 
black holes. Assessing the fraction of globular clusters and stel- 
lar systems with intermediate-mass black hole, the masses of the 
black holes and their environments, will help us to understand 
not only the dynamical processes in these stellar systems, but 
maybe also the formation of supermassive black holes. 
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